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Abstract 

In this paper, a multi-hop (MH) decode-and-forward (DF) multiple-input multiple-output (MIMO) relay network has been 
studied. To consider a more realistic scenario, Full-Duplex (FD) operation with Relay Self-Interference (RSI) is employed. 

Assuming that the MIMO channels are subject to Rayleigh fading, a simple and compact closed-form outage probability 
expression has been derived. The key assumption to derive this result is that the mutual information of each channel could be 
well approximated by a Gaussian random variable. In order to obtain the resultant outage probability, a new excellent accurate 
approximation has been obtained for the sum of Wishart distributed complex random matrices. 

Numerical Monte Carlo simulations have been performed to validate our result. These simulations have shown that, for low 
and medium interference regime, FD mode performs better than Half-Duplex (HD) mode. On the other hand, when RSI increases, 
HD mode can outperforms FD mode. 


Index Terms 

Multiple-input multiple-output, relay network. Relay Self-Interference, outage probability, Wishart matrices. 

I. Introduction 

T he demand for wireless communication technologies in the recent years seems to be unstoppable. The next promising 
history landmark will be held by 5G networks delivering massive capacity and unforeseen ways of connectivity. In 
order to achieve higher data rates, one technique that has been considered is Full-Duplex (FD) communication. FD technique 
increases the spectral efficiency for wireless systems HI since it allows devices to transmit and receive simultaneously in the 
same frequency band. However, in practice such a gain is not attainable as FD nodes suffer from Self-Interference (SI) 0 . 
therefore its performance is severely deteriorated. 

Motivation: Since FD is being considered as a promising technology in conjunction with well established techniques as 
MIMO and relay, herein, a investigation about multi-hop MIMO relay network has been conducted. In order to assess network 
performance, the outage probability metric has been used. Moreover, a comparison between FD mode and HD mode has been 
carried out to establish under what conditions a scheme is more advantageous than the other. 

Prior Related Research: The major problem to consider a network where every node is subject to Rayleigh MIMO 
channel is the mathematical complexity to manipulate random matrices. One way to reduce this complexity is to perform 
selection/combining techniques such as transmit antenna selection (TAS), space-time block codes (STBC), or maximal-ratio 
combining (MRC). Those techniques turn the problem more tractable since the final channel gains are characterize by vectors. 
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Fig. 1. System Model of multi-hop MIMO relay network, where the source (Tq), equipped with Mj^ transmitting antennas, communicates with the destination 
(Tk+i), equipped with receiving antennas, thi'ough K relays, equipped with and transmitting and receiving antennas, respectively. The solid 

lines show the desired signal, and the dashed lines show the possible RSI. 


instead of matrices a- d. Here, no selection/combining techniques were employed, therefore, the full capacity of MIMO 
system is attained. 

In im, 03, they have investigated the performance of relay network under co-channel interference and additive white 
Gaussian noise (AWGN). In ID and ifT^ . the relay is assumed to be interference-limited, thus AWGN is not considered. In 
na, if the interference levels are below certain thresholds, they treat interference as noise at the receivers. As stated in na, 
signal-to-interference-plus-noise ratio (SINR) may not be so great when the interferers also transmit signals at a power level 
similar to the source. In this case, high signal-to-noise ratio (SNR) approximation, often used in outage analysis, is not enough 
to provide an accurate outage expression. 

Contribution: Our main contribution is to derive a simple and compact closed-form expression for the outage probability. 
This expression is general for any number of antennas and SNR values. In the derivation of this outage probability, we have to 
deal with linear combination of Wishart random matrices, i.e., aWi + 6 W 2 , which is not Wishart distributed, except for the 
trivial case a = b.ln order to circunvent this problem, we propose an extreme accurate approximation based on an upper Ub and 
lower bounds Lb- Another contribution is to show that the mutual information for the MIMO relay with relay Self-interference 
(RSI) can be well approximated by a Gaussian distribution. This closed-form expression allows the performance evaluation of 
a MIMO relay network, providing a reliable way to compare FD and HD modes. To the best of our knowledge, no similar 
results can be found in the literature. 

This article is organized as follows. Section II and Section III present the system model and the outage probability expression, 
respectively. The performance evaluation is done in Section IV. Finally, Section V brings the work’s conclusions. 

H. System Model 

The scenario under consideration is described in ED, where all nodes operate in the single-frequency mode and are subject 
to relay Self-interference. Also, end-to-end communication always consists of multi-hops, i.e., direct communication between 
the source and the destination is not possible. The system diagram is partially described in ifTOl without the feedback link and 
with addition of RSI. It consists of one source (Tq), K relays (Tk), and one destination (Tk+i) as in Fig. [T] We also assume 
that only the last node, Tk+i, is not subject to RSI, since only reception occurs. The communication protocol used by any 
relay, Tk, is DF. 

Relay Self-Interference (RSI) may occurs when Tk is operating in FD mode, since it transmits and receives simultaneously. 
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The discrete-time received signal Nk x 1 at node from T^-i can be written as IflTl 


yk = + Zfc 


( 1 ) 


where x^, x^i are Mk x 1, and Mk+i x 1 transmitted signals from T^-i and nodes. The power constrains on the transmit 
signals are E l|[xj,Xfc] ^ Mk- The parameters rjk and pk, are related to the SNR ifTsl . pk = SNRk/Mk, Pk = SNRk+i/Mk+i, 


where SNRk and SNRk+i are the normalized power ratios of x^ and x^+i to the noise at each antenna of Tk node, 
respectively. Since the last node is not subject to RSI, pK = 0. The variable Zfc is independent Nk x 1 circularly symmetric 
complex Gaussian noise vector with distribution CAf{0,lNk), and uncorrelated to x^.. I^Vj, is the identity matrix of order Nk- 

The matrix Hfc is Nk x AR and represents the channel gain, as depicted in Fig. [T] Flat and spatially uncorrelated Rayleigh 
fading MIMO channels have been considered. So, is a random and independent matrix. The entries of each matrix are 
independent and identically distributed (i.i.d.) complex Gaussian variables, with zero-mean, independent real and imaginary 
parts, each with variance cr^. The matrices are known at receiver node only (receiver channel state information-CSIR). The 
matrix Hfc (/c = 1, • • • , A" — 1) is a Nk x Mk+i matrix with the same distribution as Hfc, and represents the RSI matrix. 

Define a Wishart matrix Wfc as HD 


Wfc = 


HfcH[ Nk^Mk 

HiHfc Nk > Mk 


( 2 ) 


where “f” stands for the conjugate transpose. Then Wfc is a to x m with p degrees of freedom random non-negative definite 
matrix, and thus has real, non-negative eigenvalues, with p = max(Mfc, Nk), m = min(Mfc, Nk). The distribution law of Wfc 
is called the uncorrelated central Wishart distribution, and 


Wfc- 


( 3 ) 


When Tk is in HD mode, the transmissions are orthogonal in time and therefore no RSI is generated. In this case, the 
received signal at Tk node can be written as 

Yk = v^HfcXfc -f Zfc. (4) 


III. Outage Probability 

To evaluate performance of a MH MIMO relay system, it’s needed to know if the system is able to convey a certain data 
rate, denoted here by TZ, and measured in bits/s/Hz. The mutual information (I) gives the dependence between the input and 
output of the channel, and express the amount of information that can be carried out by the channel. Consequently, the outage 
probability is the probability that the attempted T is under certain TZ, i.e., V{TZ) = Pr{I < 7^} m. 

A MIMO multiple-access channel (MAC) can be assumed since there are two transmitters and one receiver, the first 
transmission occurs from Tk-i to Tk and the second from Tk to Tk (RSI), as can be observed in ([U ifTTl . With this assumption. 


*E['] denotes the expectation operator. 
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the following rates can be achieved 11211 


^kjd ^ log 2 [det + /OfcWfc)] 

^kjd + Ikjd ^ log 2 [det (IAr^^ + rjkWk + PfeW^)] . 


(5) 

( 6 ) 


where Itjd denotes the mutual information between transmitter k and its self receiver k, and denotes the mutual 

information between transmitter k—1 and receiver k. Using the logarithm properties, and that det(AB) = det(A) det(B) and 
det(A~^) = det(A)“^, where A and B are square matrices of same dimension, and assuming that the network is operating 
in one of the corners points of the MIMO MAC capacity region mil, then 


Ikjd 


l0g2 


det 



PfcWfc \ 
Pk^k + IjVfc / 


(7) 


As a multi-hop case is being treated here, the outage probability in each of the A-hops should be evaluated. And the overall 
outage is determined by the weakest hop. Under the assumption that the hops are subject to independent fading, the outage 
probability can be written as IH 


'Pfd{TV) = 1 — PT{Iijd > TZ,... ,lK,fd > Tl} 

where Vxk /d(7^) is the cumulative distribution function (CDF) of Ik.fd- For the MIMO case, the authors in ll22ll have shown 
that I could be well approximated by a random variable with Gaussian distribution. Our assumption is that the same holds 
for the quotient of matrices in the case of considering the RSI. As will be verified in the numerical results section, this 
approximation proves to be excellent. 

For a Gaussian random variable X, the following relation is valid 1 — Vxix) = Qx{x), where Q(-) is the Gaussian 
Q-function given in ll2^ . Using this in ®, the closed-form expression for the outage probability can be written as 

K~\-l / '1D\ 

Vfd (A) = 1 - n QxFf. ( ~ ) • (9) 

k=l ^ ^kjd / 


In order to calculate the outage probability in (|9]l, it is necessary to determine Pkjd = ^[^kjd] and cr^ = E 


-Cl. 


kjd 


k-ljd- Note that, Ikjd in 0, can written as Ikjd = log 2 [det (l^v^ + -f p^W^)] - log 2 [det (lx,, +PfcWfe)]. To 

characterize this random variable, the distribution of rjk'Wk + Pk'^k would be necessary. Unfortunately, the distribution of 
linear combinations of Wishart random matrices is not known, except for the trivial case, rjk = Pk- 

The problem can be simplified using the result of ll24l that states: Let Wfe and Wfc be Hermitian m x m matrices with 
eigenvalues ai ^ a 2 ^ ^ am Pi ^ P 2 ^ Pm> respectively. In particular, if am + Pm ^ 0 (which is certainly 

true since both Wfc and 'Wk dre positive semi definite m) then 


Lb < det [ixk + PfcWfc -I- pfcWfc) < Ub 


(10) 
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where 

Lb = log 2 Y\T=1 (1 + PkOii + PkPi) 


( 11 ) 


Ub = log 2 nlLl (1 + PkCti + rjkPm+l-t) 


are the lower and upper bounds, respectively. 

A constant factor 1 has been inserted in (fTTTi to account for identity matrix present in (fTOl i. Also pk and pk have been 
inserted to account for the signal power at receiving antennas. Since Wishart matrices defined in (O obey ( fTOl i. the following 
approximation has been proposed 

log 2 [det {In, + VkWk + PfcWfc)] (12) 


Fig. in illustrates the distribution of log 2 [det {In, + + pfeW^)], considering four different combinations of rjk and 

Pk- Two important conclusions can be drawn from this figure: 1) the Gaussian approximation is perfectly acceptable; 2) the 
approximation given in (fT^ proves to be excellent since its difference from the exact distribution is almost indistinguishable. 
Therefore, the mean and variance needed are given by 


Pk,fd = E[IkJd] 

= E [log 2 [det {In, + Vk'^k + PfeW^)]] 
- E [log 2 [det (Iat^ + PfcWfe)]] 


Ri E 


Lb -\-Ub 

2 


E[A] 


(13) 


^kjd — lE[Tfc Jd] — pfejd 


E 


Lb + Ub 


-A 


2 

Pkjd 


(14) 


where 

m 

A = E [log 2 [det {In, + PfeW^)]] = ^log 2 (1 + pkoa) (15) 

as given in HD. In order to evaluate the mean E[-] with respect to the ordered eigenvalues ai and/or Pi, the joint probability 
density function (JPDF) of ordered eigenvalues ai and Pi should be used. The JPDF of ordered eigenvalues Pi ^ P 2 ^ ^ Pm 

of Wfe is given by 

f{Pi, ...,P^)= det (16) 

where K{x,y) =Y.'^i^i{^)^i{y)’ Lf{P) = Lexp{P)p-<^-^{e-^P'^+^), pi{P) = [i\/{i +L\{P)pie^ and d = p - m. 
The JPDF of ordered eigenvalues ai ^ a 2 ^ ^ Q^s of is the same as in (fTbl l. 
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n (bits/s/Hz) 

Fig. 2. Exact (marker points) and approximate (solid line) distribution of log 2 [det (Ijv^ + + PfcWj.)] for four different sets pairs of (? 7 fc, pj.). 



Fig. 3. Outage probability comparison between full-duplex an d half- duplex for the 3-hop ne twork, with SNR)^ = 20dB, SNR 3 = SNRk — 17dB, and 
Mh = = Nh = 2, for A: = 1, 2, 3. The RSI values are SNR^^ = SNR^. — 5dB and SNRi^ = RNR^; — 12dB. Markers and solid lines show the 

simulated and analytical values, respectively. 


IV. Numerical Results 

In this section, analysis for the 3-hop case, i.e., K = 2 has been done. It has been considered that the signal power associated 
with all hops is the same, and all nodes operate with the same number of antennas. Since the weakest channel determines the 
capacity, all the three channels have same parameters. 

Fig. [3 shows the outage probability for FD and HD modes with Mp = M^+i = = 2 (MIMO 2 x 2), for k = 

1,2,3, and SNRk = 20dB. For the HD mode, we have used the following expression for the mutual information Ik,hd = 
i log 2 [det (Iatj, + lyfcWfc)]. Notice that the better performance happens when there is no RSI at relays (diamond). The same 
curve is obtained when the RSI attenuation is grater than 35dB. 

Two other RSI scenarios have been considered. One with medium RSI, 12dB attenuation (square). And one with strong RSI, 
5dB attenuation (circle). Finally, HD mode scenario is also presented (triangle). It’s possible to see that HD mode outperforms 
FD mode for the last scenario (5 dB attenuation of RSI). 

More importantly, the analytical results and the simulations match perfectly, which validates our Gaussian approximation 
proposed and also the use of the lower and upper bounds. The numerical simulation was performed with 10^ channel realizations. 
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V. Conclusion 

In this paper, a simple and compact closed-form outage probability for a K relay network operating in half-duplex and 
full-duplex modes has been presented. A more realistic scenario, where the relay nodes suffer from RSI was considered. 
Our assumption is that the mutual information could be well approximated by a Gaussian random variable. Wishart matrices 
properties have been used to determine the mean and variance of the mutual information. A comparison between full-duplex 
and half-duplex was conducted. A perfect agreement between our expression and Monte Carlo simulations results has validated 
our analysis. 
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